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Abstract.  A  large  collection  of  white  dwarf  (WD)  spectra  from  Sloan  Digital  Sky  Survey 
(SDSS)  data  release  4  (DR4)  WD  catalog  data  allowed  us  to  obtain  a  luminosity  function  (LF) 
for  the  hottest  WDs.  The  LF  was  calculated  basing  on  spectra  of  the  WDs  from  a  special 
class  of  SDSS  objects  called  hot  standards  (HS),  for  which  the  WD  sample  completeness  can 
be  estimated.  WD  distances  were  determined  from  the  observed  and  absolute  SDSS  g  filter 
magnitudes  derived  from  WD  stellar  atmosphere  model  fits  to  the  SDSS  spectra.  The  result 
LF  covers  -0.5  up  to  7  absolute  bolometric  magnitude  range.  It  shows  a  slight  disagreement 
at  the  cold  end  comparing  to  the  literature  LFs  but  it  can  be  due  to  some  problems  with  HS 
WD  sample  completeness  estimations.  It  is  clearly  visible  that  the  LF  has  unexpected  plateau 
between  1.5  and  4  absolute  bolometric  magnitudes.  There  is  nothing  similar  visible  in  a  LF 
built  only  for  non-DA  WDs.  The  plateau  might  be  real  but  we  are  concern  that  its  existence 
can  be  connected  to  the  problems  with  spectra  calibration  of  the  blue  objects  in  SDSS  DR4 
data  and  lack  of  good  stellar  atmosphere  models  for  the  hottest  DA  white  dwarfs. 


A  sample  of  WDs  used  to  calculate  hot  WD  LF  was  selected  from  DR4  WD  catalog  (Eisenstein 
et  al.  2006).  The  catalog  contains  spectroscopically  identified  WDs  from  SDSS  (York  et  al. 
2000)  DR4  data  (Adelman- McCarthy  et  al.  2006),  which  covers  4783  square  degrees  of  the  sky 
with  spectroscopic  observations.  However  not  all  WDs  from  the  catalog  were  used  to  build  LF 
but  only  those  marked  as  a  special  SDSS  class  of  objects  called  HS  (Krzesinski  et  al.  2004), 
which  SDSS  uses  to  calibrate  the  observations.  Shortly,  the  HSs  are  a  subclass  of  SDSS  objects 
which  fulfil  the  following  criteria:  their  observed  g  filter  magnitude  is  greater  than  14  mag  and 
dereddened  g  magnitude  is  less  than  19  mag  (see  Schlegel  et  al.  1998).  Their  dereddened  color 
differences  ( u  —  g) o  and  ( g  —  r)o  falls  between  -1.5  and  0  magnitudes.  Objects  fulfilling  these 
requirements  are  allocated  fibers  for  spectroscopy.  In  addition  to  above  HS  criteria,  a  candidate 
for  HS  has  to  pass  some  isolation  and  quality  requirements  to  receive  a  fiber  for  spectroscopic 
observations.  These  requirements  are  very  restrictive  and  largely  decrease  HS  completeness. 
Never-less,  HSs  are  the  most  complete  of  all  stars  observed  in  SDSS  survey  because  of  their 
priorities  to  be  observed  over  all  other  star  categories  and  its  completeness  can  be  estimated. 
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The  SDSS  HS  color-color  constrains  restricting  the  sample  impose  a  WD  selection  effect  by 
cutting  out  more  DA  than  DB  class  of  WDs  from  the  sample.  In  order  to  avoid  that  effect  we 
choose  to  restrict  our  HS  WD  sample  in  the  temperature.  Therefore,  for  the  analysis  only  WDs 
with  Teff  greater  than  23500  K  and  log  g  >  7.0  were  used. 

To  calculate  our  WD  LF  a  general  inverse  maximum  volume  (l/VmSQC)  formula  introduced  by 
Schmidt  1968  was  used.  In  the  formula  I4nax  describes  a  maximum  volume  in  which  a  star  can 
contribute  to  the  LF : 

^max  =  -3-  /left  Cg  Ch  S  ( rm  a,x  —  rmin)e  250 

where,  /3  is  the  ratio  of  the  SDSS  DR4  spectral  observation  coverage  of  the  sky  (equals  4783, 
Adelman-McCarthy  et  al.  2006)  to  the  total  sky  surface  in  square  degrees.  Since  SDSS  HS  fiber 
allocation  procedure  strongly  depends  on  the  object  Galactic  latitude  b  and  object  brightness  in 
SDSS  g  filter  (Eisenstein  et  al.  2006,  section  5.6)  our  WD  sample  completeness  coefficient  will 
not  be  a  constant  value.  Therefore  two  coefficients  Cb  and  cg  defined  as  linear  functions  of  the 
object  latitude  b  ( q, )  and  its  brightness  ( cg )  were  introduced  to  take  care  of  the  incompleteness 
due  to  the  fiber  allocation.  An  additional  restriction  comes  from  HS  image  quality.  Eisenstein 
et  al.  2006  estimates  that  only  77%  of  photonretricaly  targeted  HS  objects  pass  the  stringent 
HS  photometric  quality  flags  and  the  third  completeness  coefficient  chs  =  0.77  in  the  equation 

\z\ 

above  reflects  that.  Finally,  the  e  250  factor  in  the  formula  takes  into  account  exponentially 
decreasing  WD  density  with  their  distance  z  (in  parsecs)  from  the  Galactic  plane.  The  Galactic 
scale  height  was  assumed  to  equal  250  pc  (the  same  as  in  Harris  et  al.  2006). 

The  distances  to  WDs  were  determined  from  observed  and  absolute  SDSS  g  filter  magnitudes 
derived  from  WD  stellar  atmosphere  model  fits  to  the  SDSS  spectra.  However,  one  has  to  keep 
in  mind  that  the  DR4  WD  catalog  model  fits  show  systematic  errors  which  affect  DAs  with 
effective  temperatures  greater  than  30  000  K  (Eisenstein  et  al.  2006).  The  trend  is  that  these 
model  fits  to  SDSS  spectra  tend  to  be  higher  in  temperature  than  earlier  literature  fits  and  that 
converts  into  brighter  object  in  bolonretric  magnitudes  (Mboi)-  Therefore  we  derived  a  parabolic 
temperature  correction  function  for  the  DR4  WD  catalog  Teg  and  apply  it  to  the  whole  HS  DA 
sample.  Also,  since  DR4  WD  catalog  Teg  and  log  g  fits  were  calculated  with  LTE  models,  we 
used  Napiwotzki  et  al.  1999  NLTE  corrections  to  the  log  g  and  Teg  of  all  DAs  with  Teg  > 
40  000  K.  Taking  both,  the  temperature  and  NLTE  corrections  into  account,  we  calculated  the 
LF  shown  in  figure  1. 

In  the  figure,  bolometric  magnitude  bins  are  centered  on  half  magnitudes  and  bin  widths  are 
equal  1  magnitude.  For  comparison  we  plotted  WD  LF  by  Harris  et  al.  (2006)  based  on  SDSS 
Data  Release  3  (Abazajian  et  al.  (2005)  in  the  same  figure.  As  one  can  see,  the  cold  end  of 
our  LF  is  slightly  above  the  hot  end  of  Harris  et  al.  LF  and  the  displacement  is  greater  than 
the  error  bars.  The  largest  influence  on  the  LF  position  in  such  diagrams  is  due  to  the  WD 
sample  completeness  estimation.  If  that  were  the  case,  then  our  completeness,  which  basically 
follows  the  recipe  given  by  Eisenstein  et  al.  (2006),  would  be  underestimated.  To  investigate 
the  problem  we  had  undertaken  our  own  HS  WD  sample  completeness  analysis,  independent 
on  Eisenstein  et  al.  (2006).  At  this  moment  however,  we  are  not  able  to  draw  any  certain 
conclusions  about  the  completeness  estimations  of  the  HS  sample. 

Farther,  hotter  part  of  our  LF  in  figure  1  goes  gradually  down,  then  became  horizontal 
between  1.5  and  4  magnitudes  which  is  rather  an  unexpected  behavior.  Since  the  LF  is  mostly 
populated  by  DA  WDs,  the  flat  part  of  the  LF  (or  a  plateau)  is  likely  intrinsic  to  DAs  only.  To 
make  sure  that  this  is  the  case  we  made  separate  plots  (not  shown  here)  of  DA  and  non-DA  LFs, 
which  confirmed  that  the  plateau  exists  only  in  DA  LF  and  not  in  non-DA  one.  By  separating 
DA  and  non-DA  LFs  we  could  also  demonstrate  that  there  is  a  sudden  drop  in  the  space  density 
of  non-DA  WDs  near  Mg0\  =  2  magnitudes  (which  corresponds  to  80  000  -  90  000  K)  where  the 
hottest  part  of  DA  LF  begins.  The  meaning  of  that  is  such  that  we  actually  see  the  evolutionary 
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Figure  1.  Hot  WD  LF  plot  for  the  HS  sample  and  a  bin  width  of  1  mag  centered  on  a  half 
magnitudes.  The  galactic  scale  height:  250  parsecs.  Note  that  Harris  et  al.  (2006)  LF  bin  width 
is  0.5  magnitude. 


transformation  of  non-DA  WDs  into  DA  ones.  During  the  evolution  temperatures  of  the  WDs 
are  gradually  decreasing.  As  a  consequence,  approximately  at  90  000  K,  the  effectiveness  of  the 
convective  mixing  in  their  atmospheres  decreases  allowing  for  heavier  elements  to  settle  down. 
That  way  a  pure  hydrogen  WD  atmosphere  is  created  and  a  non-DA  WD  becomes  a  DA  one. 
We  start  to  see  first  DA  WDs  in  the  DA  LF  while  at  the  same  time  the  space  density  of  non-DAs 
drops  rapidly  in  the  non-DA  LF. 

At  the  cold  end  of  the  non-DA  LF  we  do  not  see  any  unexpected  variations  but  it  is  worth 
to  notice  that  there  is  no  gap  or  any  significant  drop  of  non-DB  LF  within  DB  gap  temperature 
region.  This  is  due  to  recently  found  hot  DB  WDs  (Eisenstein  et  al.  2006)  which  filled  up  the 
gap. 

In  a  summary,  the  small  discrepancies  between  our  and  literature  LFs  can  be  easily  corrected 
by  changes  in  the  HS  sample  completeness  estimation.  From  the  other  hand  our  LF  was 
calculated  based  on  distances  determined  from  the  atmosphere  model  fits  to  the  SDSS  spectra. 
That  seems  to  be  more  accurate  approach  then  any  statistical  proper  motion  method  used  to 
built  the  LFs  in  the  literature.  Therefore  further  analysis  of  that  problem  is  required.  A  new 
light  might  bring  our  own  (independent  from  Eisenstein  et  al.  2006)  HS  sample  completeness 
estimation  from  SDSS  targeting  code  but  it  is  still  under  investigation. 

For  the  moment  we  failed  to  explain  the  reason  for  the  DA  LF  plateau  presence  with  any 
physical  cause  but  we  incline  to  say,  that  the  temperature  and  NLTE  corrections  we  applied 
should  have  removed  most  of  the  problems  with  DA  WD  temperature  and  log  g  determinations. 
In  that  case  the  plateau  we  see  in  the  LF  is  real.  However,  one  has  also  keep  in  mind  that 
there  are  other  difficulties  connected  to  the  temperature  determination  of  the  hottest  DAs  from 
the  atmosphere  model  fits.  For  example,  metals  are  not  included  in  the  NLTE  DA  atmosphere 
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models  and  their  influence  on  the  temperature  determination  from  the  models  might  be  of  order 
of  NLTE  corrections  (Barstow  et  al.  1998).  Therefore,  farther  investigation  of  the  hottest  part 
of  the  DA  LF  is  required. 
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